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Abstract 

In this paper, we construct stationary classical solutions of the shallow water equa- 
tion with vanishing Froude number Fr in the so-called lake model. To this end we 
need to study solutions to the following semilinear elliptic problem 

r_ £ 2 div (^) = 6 ( u _ g i og i)p ; mfi7 

1 u = 0, on dVL, 

for small e > 0, where p > 1, div ( ) = and Q C M 2 is a smooth bounded domain,. 



We showed that if has m strictly local minimum(maximum) points z^, i = 
1, ■ ■ ■ ,m, then there is a stationary classical solution approximating stationary m 
points vortex solution of shallow water equations with vorticity YllLi "fjjy^- Moreover, 
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strictly local minimum points of on the boundary can also give vortex solutions for 
the shallow water equation. As a further study we construct vortex pair solutions as 
well. 

Existence and asymptotic behavior of single point non-vanishing vortex solutions 
were studied by S. De Valeriola and J. Van Schaftingen in [9]. 

Mathematics Subject Classification: 35J20, 35J40, 35J60, 35R45. 
Keywords: Shallow water equation; Free boundary; Stream function; Vortex. 

1 Introduction and Main Results 

We consider fluid contained in a basin by a uniform gravitational acceleration g and fixed 
vertical lateral boundaries(i.e. no sloping beaches). Suppose that (x, y) is horizontal spatial 
coordinate which is confined to a fixed bounded domain Q with boundary dQ.The vertical 
coordinate is chosen so that the mean height of the fluid's free upper surface is at z = 0. Let 
z = —b(x, y) give the fixed bottom topography, so b is a strict positive function over Q. Let 
z = h(x, y) be the free upper surface. We assume that both b and dQ vary over distances L 
which are large compared to the mean depth B, that is, the ratio 5 = j- is small. 
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Let u and w denote the horizontal and vertical components respectively of the fluid veloc- 
ity. We will consider only those motion for which u,w and h each vary in (x, y) over distances 
L, in other words, we will make the long- wave approximation. The "Froude number" is de- 
noted as Fr = -r=s, where U is the characteristic magnitude of u. We will consider the case 
of small "Froude number" Fr and h is small compared to B. In such cases, from [U El SHU] 5 
the leading-order evolution of u(x, y, t) and h(x, y, t) will be governed by equations that have 
the non-dimensional form 

'<9 t v+(v V)v = -V/i, 
div(fov) = 0, 

where V is the horizontal gradient. Since these equations apply to a domain which is shallow 
compared to its width and whose free surface exhibits negligible surface motion, they are 
called the 'lake' equations(see [3], for instance). 

The first equation in (11. ip can be rewritten in terms of the vorticity u = V x v as 

<9(V + u x v = —V ( + h 



This model is analogous to the two-dimensional Euler equation for an idea incompressible 
fluid and has been recently studied by many authors. For instance, see |TJ [3j HJ [15] and the 
references cited therein. 

Recently, De Valeriola and Van Schaftingen [9] studied the desingularization of vortices 
for ( 11. ip with stream function method, which consists in observing that if i\) satisfies 

-div {^j = bfW) 

for / G C\R), then v = ^ and h = -pty) - ^ with F(s) = J Q S f(s)ds form a stationary 
solution to the shallow water equation. Moreover, the velocity v is irrotational on the set 
where fiif)) = 0. It is easy to see that if ipQ satisfies div (^y 2 -) = 0, then v = cur ^° is an 
irrotaional stationary solution of ( ll.ip . In [9], they studied the asymptotics of solutions of 



-£ 2 div(^) = bifj p + , intl, 
ip — ipo In - , on dQ 



where p > 1, Q C M 2 is smooth bounded domain. 

To obtain their results, De Valeriola and Van Schaftingen investigated the following 
problem 

= b(u - q e ) p + , in fi, 




on oil, 

where p > 1, q = —ipo, q e = q In - , Q C M 2 is smooth bounded domain. 

More precisely, they first obtained the existence of solutions by using mountain pass 
lemma and studied the asymptotic behavior of solutions by giving exact estimates to the 
upper and lower energy bounds of the least energy solutions. As a consequence, they obtained 
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that the "vortex core" shrinks to a point Xq which is the minimum point of y. However, it 
is hard to apply their method to construct multiple vortices for fll.ip . 

Motivated by [5], our goal in this paper is to construct multiple stationary vortices for 
shallow water equations. More specifically, we want to find some high energy solutions whose 
"vortex core" consists of multiple components which shrink to several distinct points in fl 
as e — > under some additional assumptions on ^P-. 

Our main results in this paper can be stated as follows: 

Theorem 1.1. Suppose that fl C M 2 is a smooth bounded domain. Let b G C 1 (f2) ; -0 O G 
C 2 (fl) be such that div(^p) = and let Vq = curh/V //info 6 > and sup Q < 0, then for 

any given strictly local minimum (maximum) points Z\, ■ ■ • ,z m of^P-, there exists Eq > 0, such 
that for each e G (0, Sq), there exists a family of solutions v £ G C^f^M 2 ) and h e G C l (fl) of 

div(bt; e ) = 0, in fl, 

[v £ ■ V)v E = -Vh E , infl, 
v £ ■ n = v ■ nln ~, on dfl, 

where n is the unit outward normal. Furthermore the corresponding vorticity uj £ := curl v £ 
satisfying 

suppu £ C U™ =1 B(z ij£ , Ce) forz ij£ G fl, i = 1, • • • ,m 



and as e — > 0, 

2-K1p Q (Zi 



in 



E 



1b 2 

The next result shows that strictly local minimum points of -P- on the boundary dfl can 
also give vortex solutions for f ll.ip . 

Theorem 1.2. Suppose that fl C M 2 is a smooth bounded domain. Let b G C 1 (n) ; -0o G 
C 2 (n) be such that div(^g^) = and let v = curl^ . //inf^fe > and sup^^o < 0, then, 

for any given strictly local minimum points Z\, - • • ,z n of^P- on the boundary dfl, there exists 
6q > 0, such that for each e G (0,£o); there exists a family of solutions v £ G C 1 (n,IR 2 ) and 

K g c\fi) of 

div(bt; e ) = 0, in fl, 

(v £ ■ V)v £ = -Vh £ , in CI, 

v £ ■ n = Vq ■ nln i, on dfl, 

where n is the unit outward normal. The corresponding vorticity u £ := curl v £ satisfying 
suppu £ C \Ji =1 B(z i}E , Ce) forz i)£ eft, i = 1, • • • , n 

and as e — > 0, 

E24(z«) 



£ 

n 



^ Kb) 
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Moreover, 

A . „ /In I In el \ 1 

\Zi,e ~ Zi\ < C I I , dM^^e, dSZj > | ln£ | Q > 

where a is a positive constant. 

Remark 1.3. If 'fy- has strictly local minimum points in Vl and on the boundary dVl, then 
there is, from Theorem \l.l\ and \l.S\ a stationary solution of the shallow water equation such 
that its vorticity set shrinks to corresponding strictly local minimum points . 

It is worthwhile to pointing out that although the structure of shallow water equations is 
very analogous to that of two dimensional Euler equations for an ideal incompressible fluid, 
the position of vortex for ( II. ip exhibits a striking difference with that of the Euler equations. 
The position of vortex for Euler equation is closely related to Kirchhoff-Routh function. The 
interested reader can refer to [3 El [161 EH] f° r more results on this problem. 

Theorem 11.11 and Theorem 11.21 are proved via the following results concerning problem 

Theorem 1.4. Suppose that Q C M. 2 is a smooth bounded domain. Let b G C 1 (f2), q G 
C 2 (Cl), infnb > and inf^ q > 0. Then, for any given strictly local minimum (maximum) 
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points Zi, ■ ■ • ,z m of ^, there exists Eq > 0, such that for each e G (0, Eo), ( 11 .3j) has a solution 
u £ , such that the set fl £ = {x : u £ — gin - > 0} has exactly m components fl £t i, i = 1, • • • , m 
and as e — > 0, each Q £t i shrinks to the point Z{. 

Theorem 1.5. Suppose that Q C M 2 is a smooth bounded domain. Let b G C 1 (fi), q G 
C 2 (fl), iniftb > and inf^ q > 0. Then, for any given strictly local minimum points 

2 

, z n of 2- on the boundary d£l, there exists e > 0, such that for each e G (0,e ); 
(II. 3p has a solution u £ , such that the set Q £ = {x : u £ — gin | > 0} has exactly m components 
f2 £) i, i = 1, ■ ■ ■ ,m and as e — ¥ 0, each fl £ j shrinks to the point z^. 

Not as in [9] where (11.31) is investigated directly, we prove Theorem 11.41 and Theorem 

11.51 by considering the following equivalent problem of ( 11.31) instead. Set 5 = e (in ^) 2 , 
w = jif^r, then (jl.3p becomes 

(S 2 div(^)=b(w-qr + , inO, 
1 u — 0, on dVt. 

We will use a reduction argument to prove Theorem 11.41 and Theorem 11.51 To this end, 
we need to construct an approximate solution for (jl.4p . For the problem studied in this 
paper, the corresponding "limit" problem in M 2 has no bounded nontrivial solution. So, 
we will follow the method in [8] to construct an approximate solution. Since there are 
two parameters 5, e and b in problem ( 11.41) . which causes some difficulty, we must take 
this influence into careful consideration and give delicate estimates in order to perform the 
reduction argument. 

We will also apply the above idea and techniques to construct vortex pairs to shallow 
water equations in section 5, which has never been addressed before. 
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As a final remark, our results seem connected with the work of Wei, Ye and Zhou [T7J 
[TBI US] on the anisotropic Emden- Fowler equation 

div(a(x)Vw) + e 2 a(x)e u = 0, in f2, 
u = 0, on dQ. 

They constructed (boundary) bubbling solutions showing a striking difference with the isotropic 
case(a = constant). Moreover, we point out that problem (11.41) can be considered as a free 
boundary problem. Similar problems have been studied extensively. The reader can refer to 
[3 El [7J HI [8j [TTJ [13j HH] f° r more results on this kind of problems. 

This paper is organized as follows. In Section 2, we construct the approximate solution 
for (11. 4p . We will carry out a reduction argument in Section 3 and prove the main results in 
Section 4. In Section 5, we give some further results on vortex pairs for the shallow water 
equations. Some basic estimates that used in sections 4 and 5 will be given in Section 6. 



2 Approximate Solutions 

In the section, we will construct approximate solutions for (11.41) . 

Let R > be a large constant such that for any x G Q, Q CC B R (x). Consider the 
following Dirichlet problem: 

' -5 2 Aw = (w - af + , mB R (0), 
w = 0, ondB R (0), 

where a > is a constant. Then, (12. ip has a unique solution Ws >a , which can be written as 

w^)4 a :t s '^ ] ' N - si ' (2.2, 

[«]^r, ss < \x\ < R, 

where 4>(x) = (f>(\x\) is the unique solution of 

-A0 = ^, 0>O, ^6^(Bi(0)) (2.3) 

and ss G (0, R) satisfies 



In | 



which implies 

£-1 

2 



S\\n5\ 



— -> ( — : — ) > 0, as 6 -> 0. 
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(1 



Moreover, by Pohozaev identity, we can get that 



p+ l= 7r(p + l) ; a ^ , ( u> = 2n \ ( i/(i) 



Bi(0) ^ ^Bi(O) 
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For given b > and q > 0, let Vg^Jx) be the solution of the following Dirichlet problem 

-5 2 Av = b 2 (v-q) p + , mB R (0), 
v = 0, ondB R (0). 

By scaling, from ( 12.1ft and (12. 2ft . we obtain 




v~( x ) = r— (x) = r .. V"/ YK « J ' 1 1 - (2.5) 

S(5 < N < a. 



For any z£ fl, define V^g^ 2 (x) = V^g -(x — 2). Because V^g - does not vanish on dfl, we 
need to make a projection. Let PV 8 ^~ Z be the solution of 



-8 2 Av = P{V sM>z - q) p + , in O, 
v = 0, on <9fi, 



(2.6) 



and h(x,z) be the solution of 



-A/i = 0, in fi, 

ft, = =- In ] r . on <9fi. 

2-7T |as— «| ' 



Then 

In 



= - 7Tn9(x, z), (2.7) 



where #(x, z) = In i? + 2irh(x, z). 

We will construct solutions for (II ,4p of the following form 

m 

where Zj G Q for j = 1, • • • , m, oog is a perturbation term. To obtain a good estimate for us, 
we need to choose qs,j properly. 

Denote Z = (z\, - ■ ■ , z m ) G R 2m . In this paper, we always assume that zj G Q satisfies 

\zi — Zj\ > q l , dist(zj, d£l) > q > 0, or 

1 (2.8) 
\zj -Zj\ < rj, dist(sj,Sft) > , i,j = l,--- ,m, 1 ^ j 

where q, r/ > is a fixed small constant and L, a > is a fixed large constant. 

Let bj = b(zj) and qsj(Z), j — 1, • • ■ , m be the solution of the following problem: 

ft = + ^li9{zi, Zi) - riiG(Zi, Zj), (2.9) 
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where G(x, Zj) = In , ^ z , — (7(0;, z,). It is not difficult to see that since In j — > oo as e — > , 
( 12. 9 p is a linear system with coefficient matrix, which is a small perturbation of a positively 
definite diagonal matrix for small e. Thus we can obtain the solution (qs,i(Z), ■ ■ ■ , q$ jm (Z)) 
to ( 12. 91) . Moreover, we have 

- 1 g (*,*) ■ 

1 lnf 

For simplicity, for given Z — (z\, ■ ■ ■ , z m ), in this paper, we will use q$,i instead of q$^(Z). 
Define 

m 

Vs^-PV^^. Vs, z = J2Vs,z, r (2.10) 

i=i 

Let be the solution of 

A 2 

r _? - — H~ l <l8,i 

m s 

then we have 

1 1 /lnllnelA 

+ O ' 1 



ss,i £ 



V l ln£ l 2 / ' 



Thus, we find that for x € BL SSi { z i)i where L > is any fixed constant, 



SS,i 

s <5,i V 



R 



e 



and for j ^ i and x G B Ls&x (zi) 



So, by using (12. 9p . we obtain 
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V s>z (x) - q(x) = V St i i4s . iZ .{x) - q 5 ,i + O y j^^ 9(zi, Zi^j , x G B Lss ,(zi). (2.11) 



We end this section by giving the following formula which can be obtained by direct 
computation and will be used in the next two sections. 



dz i)h 



= < 



Q5,i %i,h Xh 

In — \x — Zj\ 2 

SS.i 1 1 



+ 0(1), 



xen\B a (zi). 



(2.12) 



3 The Reduction 



Let Vs t z be given as in (12.101) . we are to find solutions of the form V^ z + u>s,z, where u^z is 
a small perturbation(obtained in Proposition 13 .3 j) . We will show that for any given Z, there 
exists us, z such that wg > z = Vs >z + w<s, z satisfies 



5 2 

—Vw &)Z ^v - b(w s ,z ~ q)+v 
b 



0, for any v G £Tj}(fi) D W 2 ' P {VL) \ H* 



(3.1) 



where H* is a finite dimensional subspace of Hq(Q) D W 2,p (Q). In the next section, we will 
choose Z properly so that Vs t z + ^8,z is a solution of (11.41) . 

To show (13. ip . we need to study the kernel of Cw := — <5 2 div(^) —pb(Vg t z — q)+~ w- To 
do this first we need to understand the kernel of the linearized equation of 



- Aw 



w 



in 



(3.2) 



Let 



w(x) 



<K\*\), 

b'(l)\n\x\ 



M < i, 

\x\ > 1, 



where 4> is the solution of (12. 3p . then w G C^IR 2 ) is the unique solution of (13. 2p . Since 
0'(1) < and In |x| is harmonic for \x\ > 1. Moreover, since w + is Lip-continuous, by the 
Schauder estimate, w G C 2,a for any a G (0, 1). 

The linearized equation of (13. 2p at w is as follows 



— Av — piu^ v = 0, v G L c 



(3.3) 



It is easy to see that i = 1,2, is a solution of (13. 3p . Moreover, from Dancer and Yan 
[8], we know that w is also non-degenerate, in the sense that the kernel of the operator 
Lv := -Av-pw^\ v G D 1 ' 2 (R 2 ) is spanned by {JJ, ||}. 
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Let Vg t z,j be the function defined in (12.10p . Set 

F s ,z= ^u:ueL p (n), IJ^± u = U,j = l i ... , m , h = \,2\, 



and 



E s ,z = : «£ W 2 > P (Q) n ff^fi), jf A (j^J « = 0, J = 1, • • • ,m, h = 1, 2} . 
Define, for any it G L P (Q), Q^it as follows: 

i=i /i=i ^ j ' h ' 

where the constants Cj t h{j = 1, • • • , m, h = 1, 2) are chosen to satisfy 




Since J* n g^'.^ Q<5^ = 0, the operator Qs can be regarded as a projection from L P {VL) 
to F^. In order to show the existence of Cj^ satisfying (13.51) . we just need the following 
estimate ( by (I2TT2]) ): 



y n Mw* - m + i — ^ J ^- + InnTpJ (3 - 6) 



x d n 



i lne|P+ 1 Vl ln£ l p 

where d > is a constant, 5 ij7l ^ = 1, if i — j and h — h; otherwise, = 0. 
Define 



L s u = -5 2 div C-j^j ~ P b ( V s,z - <?) + 1 u - 



For the operator QsL$ we have the following lemma. 

Lemma 3.1. There are constants p > and 5 > 0, such that for any 5 G (0,<5 ] ; Z 
satisfying (I2.8p . it G with Q5L5U = o(l) m L p (fi \ U™ =1 5l S4j (^)) /or some L > /arge, 
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Proof. We will use || • || p , || • to denote || • \\lp(Q) an d || • ||i°°(n) respectively. We argue 
by contradiction. Suppose that there are 8 n — > 0, Z n = (z% !n , ■ ■ • , z m ^ n ) satisfying ( 12. 8 p and 
u n e E Sn! z n with Q Sn L Sn u n = o(l) in L p (Vt \ ^f = iB LsnJ {z^ n )), such that 



\\Q 5n L Sn u n \\ p < 



1 



n I ln£„| p 



and 1 1 w„ | |oo — 1, where and in the sequel we set s n j = ss„j to simplify notation. 

Firstly, we estimate Cj^n corresponding to u n in (13. 4p . By definition Cj^ n satisfies: 



Qs n L Sn u n = L Sn u n - 5^ c ^' n \ ~^ dV Q z Zn ) 

7=1 fe=l ^ ' 



we obtain 



For each fixed i, multiplying ( 13. 7p by 9V ^"' Z "'' , noting that 

i.h 



{Qs n L Sn u n )— = 0, 

oz i,h 



(3.7) 



dV* 



S n ,Z n ,i 



n \ dz ih 

m 2 

~- a '- h 



i L 5 n U n ) 



j=l h=l 

Using ( 12. lip and Lemma l6.1[ we obtain 



-51 A 



dVs n ,z n ,i 
n dzih 

9V Sn ,z n ,j\ dV Sn ,z n ,i 



dZi 



it 



u n L 5 



dz, 



i.ti 



U,. 



h dZi 



Ph {V 5n ,z n - q) 



Ur. 



1\ _dV* 



+ P U r . 

In 



--0 
--0 



(hi - b) (Vk,z, - q) 



dz 



i,h 

i 9V Sn> z n , 



dz 



e n In 2 | 


lne n | 




lne„ 


|p+i 



ln^lf" 1 



p-1 ^ 



dz 



i.ti 



\ne n \P 



e n In 2 | 


lne n | 




lne n 


|p+i 



Using ( I3.6p . we find that 



Ci,ft,n = O (^nln 2 | lne n |) . 
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Therefore, 



m 2 



i /dV s 



dVs n ,z n ,j 

j=l h=l ^ h 
m 2 

= 5E 5E P b2 j,n C 3,h,n (^S,-, „,Zj, „ ~ <?j,n y \ 
j=l ft=l V 

+ X 5E 2 h,n c j,h,n-g^ (y8nhn,qj,n,Zj,n ~ ^>). 



<5n }bj,niQj,n i z j 



3=1 h=l 
O 



EE 

\j=l /i=l 

2 



til 


Cj,h,n | 




In e n 


P 





In 2 


lne n | 




\ne n 


p 



in L p (fi). 



Thus, we obtain 

-^<5 n M n = Q Sn L Sn u n + O 
For any fixed i, define 



F p 


In 2 


lne n | 




lne n 


V 



o 



1 F p 



n I lneJ p 1 



«i,n(y) = Un(s nii y + 2j jn ). 



Let 



L n M = — div 



Then 



Noting that 



we find that 



2 ^ _ 

p u n r ~ 



I Lx U 



S n u n\\p- 



=0 



8 n 



1 



\nen\P- 1 

2 



As a result, 



£„Ui,„ = o(l), inL p (fi„), 
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where Q n = {y : s n4 y + z ijU eQ}. 

Since |]€fc» n,|]oo — 1> by the regularity theory of elliptic equations, we may assume that 



Ui, n -)■ Ui, in C} oc (M. 2 ) 



It is easy to see that 



s 2 . 

-jf (V Sn ,Z n (s n ,iy + Zi, n ) - q{s n ,iy + Zi,n)Y + 1 
n 

c 2 / / 1 2 i i _ I \ \ P- 1 



In 2 


lne n | 




lne n | 


2 



1 



p-1 

/ 111 Mil " I \ \ 

<P(y) + o 



In 2 


lne n | 




lne n 


2 



Then, by Lemma 16.11 we find that Ui satisfies 

— Aui — pw^~ l Ui = 0. 
Now from the Proposition 3.1 in [8], we have 



dxi dx 2 

Since 

A — )u n = 0, 

we find that 

"-^. = 0, 



dw dw 

Ui = ci- Vc 2 ^—. (3.8) 



which, together with (I3.8p . gives = 0. Thus, 

Ui, n -> 0, in C^S^O)), 

for any L > 0, which implies that -u n = o(l) on dBL Sni (zi, n ). 
By assumption, 

Q 5 „L 5nMn = o(l), in Z7(fi \ ut XJ B iSn .(^ >n )). 
On the other hand, by Lemma 16. 1\ we have 

(V Sn ,z n ~ q) + = 0, x G n \ U^B^ (*,„). 

Thus, we find that 

-div = o(l), in n \ U™^.^). 
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However, u n = on dQ and u n = o(l) on dB LSn .(z i n ), i = 1, • • • , m. So we have 

u n = o(l). 

This is a contradiction. 

□ 

Proposition 3.2. Q^L^u is one to one and onto from E$ z to F&z- 

Proof. Suppose that Q5L511 = 0. Then, by Lemma l3TTj u = 0. Thus, Q5L5 is one to one. 
Next, we prove that QsLs is an onto map from E^z to F$ z . 
Denote 

E= \u:ueH^Q), / S/^^-Vu = 0, j = I,--- ,m,h = 1,2). 

Note that % z = £ n W 2 ' p (fi). 

For any /t G by the Riesz representation theorem, there is a unique u G ^(fi), such 
that 

5 2 [ VuVy? = [ Jup, V</?G#o(fi). (3.9) 
Jn Jq 

On the other hand, from h G F$,z, we find that u G .E. Moreover, by the //-estimate, we 
deduce that u G K /2,p (fi). As a result, -u G Eg jZ . Thus, we see that Qs(— S 2 A) = —5 2 A is an 
one to one and onto map from E$,z to F$ t z- On the other hand, QsL$u = h is equivalent to 



u = <T 2 (-QfA) -1 [Tu + bh], u G E s>z (3.10) 

where 



Tu = b5 2 v\vu + pb 2 (Vs, z - q) p - 1 « + E E ("^ • 

j=i h=i ^ ^ ' 

It is easy to check that 5~ 2 (— Q$A)~ l Tu is a compact operator in E$ t z- By the Fredholm 
alternative, (I3.10p is solvable if and only if 

u = 5- 2 (-Q s A)- l Tu 

has only trivial solution, which is true since QsLs is a one to one map. □ 
Now consider the equation 

Q s L s u} = Qsh + QsRs(u), (3.11) 

where 

m b 2 f 1 \ 

h = b (Vs, z - qf + -Ef - fo)' + 52 ( V ^ V M ( 3 - 12 ) 
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and 



Rs(u) =b (V s ,z + co-q) p + -b (V s ,z - qf + - pb (V s>z ~ qf^ w. 
Using Proposition 13. 2\ we can rewrite (13. lip as 

u = G & u =: {QsL s )~ 1 Qs(ls + Rs(u)). 



(3.13) 



(3.14) 



The next proposition enables us to reduce the problem of finding a solution for (ll.4p to 
a finite dimensional one. 

Proposition 3.3. There is a 5q > 0, such that for any 5 G (0, 8q] and Z satisfying ( 12.81) . 
(13. lip has a unique solution us G Eg jZ , with 



o 



In 


lne| 


lne 


2 



Proof. It follows from Lemma 16.11 that if L is large enough, 5 is small then 

(Vs,z-q) + = 0, xen\U™ =1 B Lss .( Zj ). 



Let 



M = J E? 4iZ n{H| 00 <— 



where # > Ois a small constant. 

Then M is complete under L°° norm and is a map from to Es,z- We will show 
that Gs is a contraction map from M to M by two steps in the following. 

Step 1. Gs is a map from M to M. 

For any G M, similar to Lemma [6.11 it is easy to prove that for large L > 0, 5 small 



(V s ,z + u-q) + = 0, in O \ Uf =1 B Lss .{ Zj ). 
Note also that for any u G 

Q 5 u = u inn\Uf =1 B Lss .{ Zj ). 
Direct computations yield that 

r-2 



(3.15) 



l<p<2, 



Therefore, using Lemma [6. 1[ (13.121) and (I3.13p . we find that for any u G M, 



Qsh + QsRs(u) = o(l), in L?(Q \ U™B Lss .{ Z j)). 
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So, we can apply Lemma [3.11 to obtain 

\\{Q 5 L 5 )- l (Q 5 l s + Q s R s (oj)) Hoc < Ce-^\\ne\ p - l \\Q 5 l s + Q 5 R 5 {u)\\ p . 
Thus, for any u G M, we have 

l|G' 5 H|| 00 =||(g (5 L (5 )- 1 (g 5 / 5 + Q (5J R (5 H)|| 00 



(3.16) 



<Ce"\\ne\ p - 1 \\Q s k + Q 5 Rs{u, )Hi ,. 
It follows from f)3.5p -f l3TB|) that the constant Cj^, corresponding to u G L°°(Q), satisfies 

\dV s ,z,i 



Ic^l^Clln^r 1 ^ f 

i, h 



\u\. 



Hence, we find that the constant Cj^, corresponding to Is + Rs(u) satisfies 



\c hh \ ^Clln^r 1 ^ / 

i, h 



8V* 



6,Zj 



where 



As a result, 



\k + Rs{oj 
dV s , z 



dzi~ h 



lit n 

<ciin e rj:E/ 

i,h 3=1 JBls S,j(zj) 

<Ce l ' 1 v\ \ne\ p \\k + R 5 (uj)\\ p + Ce 2 \ lne 



\k + R s {u)\ + Ce 2 \\ne\ 



b 2 / 

is = b(Vs, z - q y + -"£i {v 6 

3=1 



S,bj,qg j,Zj 9S,j 



where 



\\Q s {I s + R 6 {lo))\\ p 
<\\1 S + R s (uj)\\ p + Cj2\ c j,h\ 



3,h 



dz j>h J 



^||I 5 || p + C||^(a;)|| 3 



1, 



Ce p 
lne IP" 1 



--C\\ls\\ p + C\\R s (u)\\ p + R(p,e) 

Ce 2 



R(p,e) 



\ne\P 



, Kp<2, 



On the other hand, from Lemma [6.11 and (12. lip , we can deduce 
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S\\p 



< 



m U2 



3=1 



C Inline! 



<C 



1 llneP 



£p In | hie | 



p-1 



, lne|p +1 

For the estimate of ||i?5(o;)||p, we have 

||^M|| P =\\b(V 5 ,z + co- qf + -b(V 6 ,z - q) P + -pb(V s , z - qf' 1 

<cMl\\(Vs, z - q f-% 



<c £p \M\oo 
~ |lne|P- 2 ' 



Thus, we obtain 



||G*HI|oo <Cs'l\ lner 1 (||Z,|| p + H^HIIp + i2(p, e)) 

! 'lne|P +1 |lne|f- 2 



< 



lne| 



2-6 



Thus, (7,5 is a map from M to M. 

Step 2. G,5 is a contraction map. 

In fact, for any Ui G M, i — 1, 2, we have 



G^i - = (Qj-La) ^ (^(^i) - ^(u^)). 



Noting that 



RsM = R sM = 0, in Q \ Uj =l B Lss .{ Zj ), 
we can deduce as in Step 1 that 



\\Gguj, - GW2IL <Ce~\ Ine^WRsM - R 5 (u 2 )\\ p 



lne:|P- 2 |lne|P" 
C , , 1„ 

-\\ n£ \l-d^ Ul _ ^2||oo < g 11^1 -U2\\oo- 
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Combining Step 1 and Step 2, we have proved that G$ is a contraction map from M to 
M. As a consequence, there is a unique oog G M such that ujs = Gs^s- Moreover, it follows 
from (15351) that 



k^5 oo < C 



In 


lne: 


In e 


2 



□ 



4 Proof of Main Results 

In this section, we will choose Z properly so that V$ y z + ^s, where us is the map obtained in 
Proposition I3.3[ is a solution of (11. 4p . 
Define 

and 

K(Z) = I(Vs,z + tu s ). (4.1) 

It is well known that if Z is a critical point of K(Z), then Vs.z + ojg is a solution of ( 11.411 . 

In the following, we will prove that K[Z) has a critical point. To do this let us first show 
that I iVs,z) is the leading term in K(Z). 

Lemma 4.1. We have 

/^W(^) + («). 

Proof. Recall that 

m 

i=i 

By the definition of K(Z) 

K{Z) =1 (Vs, z ) + 6 2 [ hv s , z Vco 5 + t ! 1M1 

P + 1 Jn 

Using Proposition 13.31 and (I3.15p . we have 

/ b [(V s ,z + oos- q)l +1 ~ (Vs, z - q)l +1 ] 
Jn 

4j\H\oo \ = f s 2 ln|lne| \ 
|lne|f J V |lne| p+2 / ' 
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On the other hand, 



52 L l vv «> zVus =%S3 fc**** - us - 52 L v l vvs ^ s 



o 



e 2 In I lne| 
|lne|P+ 2 



Finally, we estimate S 2 f n ^\Vus\ 2 - 
Note that 



™ U2 . p 



(1 \ b 2 

-vus) =b(v s , z + uj s -q)i-j2j {y& 

1 m 2 , 

+ 6 2 V-VV 5 ,z + c J> h ( 

7=1 h=i ^ 



3 

m 2 / . ., - 

<9K» 



dz 



j,h 



hence 



5 2 !\\Vu 5 \ 2 = ! b(V s , z + u 5 - q y + u; s -J2 I i(v s , 



2 



ft 



0J S 



' e 2 In 3 I In el \ ^/e 2 ln|lne| 
=0 -r: ^r 1 +0 



--0 



lne|P+ 3 y V |lne|P+ 2 
e 2 ln|lne| \ 
|lne|f+ 2 ) ' 



□ 



Proof of Theorem l.J±. By Proposition 16. 2\ we have 

V" 7152 g 2 (%) 1 ^ ^ 2 ln|lne| ^ 

K{z) -^K^ +0 h^r)- 

Since zj, • • • , z m is strictly local minimum(maximum) points of for 5 > small enough, 
there exists a neighborhood £) j of Zj, ? = 1, • • • , m, such that the reduced function K(Z) 
admits at least one critical point in O e ^. Hence, we get a solution ws for (11.41) . Let u e = 
I lnelw^ and S — e| lne| 2 ( it is not difficult to check that u £ has all the properties listed in 
Theorem 11.41 and thus the proof of Theorem 11.41 is complete. □ 

Proof of Theorem 11.51 Define 



M 



jz = (z u ■ ■ ■ , z n ) E n n : \ Zj - %| < rj, distfo, 90) G (y^p } 
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where T\ and r 2 will be determined later. 

Consider the following minimizing problem 



min K(Z). 



zeM 



There exists a minimizer Z £ for K(Z) in M.. Now, as in Theorem 11.41 we just need to verify 
that Z £ is an interior point of M. and hence is a critical point of K(Z). 
By Proposition I6.2[ we have 



.r- 

Let Z e — (z E i, • • • , z £ n ) £ M. be such that 



5 2 


In 


lne| 




lne 


2 



then 



b(z Ed ) b(zj) 
As a result, 



|z £J - %| = dist(z E|i , = j^p, 



O frj~p) ' 9{2ej,Zej) = 0(ln|lne|), j = 1, 



R7 7^ V" wS2 3 2 (%) i i n ^ 2ln l ln£ l ^ 



Note that 
we find 

and 



-lnf b(zj) 

K(Z £ ) < K(Z £ 



Q 2 ( z e,j) _ g 2 (%) < c /ln|lne| 



b (Ze,j) &(%) V l ln ^ 



g(Ze,j,z e ,j) < Chi I lne|, j = 1, ■ • • ,n, 

where C is independent of t\ and r 2 . 
Hence, for j = 1, • • • , n, we have 

|z e ,;j - %| < C f "TT~T~ ) ' dist(z £ j,dQ) > 



lne| / | lne| c 

Thus, Z £ is an interior point of M. if we choose Ti to be sufficiently large and r 2 sufficiently 
small in the definition of Ai . □ 

Proof of Theorem \l.l\ and \l.S\ . By Theorem 11.41 and II .5[ we obtain that u £ is a solution of 

(OD- 

Define for iGO, 

_ curlp £ - q £ ) 

k _ b(u £ -q £ y + +1 |v £ | 2 
(p + l)e 2 2 
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Then, v e is a stationary solution of (11. ip with curlv e = -^{u e — q e )+- 
What remains to do is just to verify, as e — ¥ 0, that 



„ m 

/ curlv £ — > 



2nq(zj) 

By direct calculations, we can obtain for e small that 
b 



curlv £ = / j^{u £ - q £ Y + 



a Jn £ 

lnd p 



f b(V s , z + u s -q) p + 

/ bdV 5 , z + Us-q)l + jr 1 -^ [ { b-k){V s , z + u s -q)\ 

i=l 6 JB LHti [zi) j=1 £ JB LsSa ( Zl ) 



^ \lne\P f t / \p „ /ln 2 |lne|\ 



In 


2 


lne| 




ln£ 





^ 27rgg,j | lne| / ln 2 |hi£| \ ^ 27rg(^ 



Therefore, the result follows. □ 



5 Further Results 



In this section, we will use the idea and techniques in the previous sections to construct 
vortex pairs for the shallow water equations. For this purpose, instead of (jl.3p . similar to 
[9], we now consider the following boundary value problem: 



-eMiv(^) = b(u - q £ f + - b(-u - q £ f + , in ft, 
u = 0, on <9ft, 

where p > 1, q = —ifio, q s = <?hi -, ft C R 2 is smooth bounded domain. 

i-p 

Set 5 = e (in -) 2 , w = r^r, then (15. ip becomes 

'-5 2 div(^) = &(w - q) p + - b{-w - q) p + , in ft, 
w = 0, 



on <9ft. 



Let 6f = 6(2^) and qf i be the solutions of the following system: 



e 



k^i £ 



1=1 £ 



£ k^j £ 1=1 £ 



(5.1) 



(5.2) 



(5.3) 
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where % = 1, • • • , m, j = 1, 
Set 

and 



Let sf { be the solution of 



op- 1 s p- 1 <p (1) = 



Then, we have 



In -§- In - 



o 



In 


lne 


lne 


2 



Thus, as in ( 12. lip , we find, for z = 1, • • • , m and j = 1, ■ ■ ■ ,n that 



and 



Similar to Proposition 16. 2\ we have the following energy expansion: 

J(V+ _ y- ) _ A , ^ + ) 



-lnf &(*+) ^ In? 
A ttJ 2 g 2 (z~) ( g(zj,z~) \ 
~t ln 7 h ( z j) \ lnf ) 



,'5 2 ln|lne|^ 



j 

From ( 15.41) . we can deduce the following result: 



(5.4) 



hid 



Theorem 5.1. Suppose that Q <Z M 2 is a smooth bounded domain. Let b G C 1 (fi) ; -0o € 
C 2 (f2) be such that div(^g^) = and let Vq = curh/V U^nb > and sup^^o < 0, then 
for any given strictly local minimum (maximum) points zf , ■ ■ • , z+, z± , ■ ■ ■ ,z~ of ^P-, there 
exists Eq > 0, such that for each e G (0,£o)? we can find a family solutions v £ G C 1 (f2,R 2 ) 
and h £ G C l (Q) of 

div(bv £ ) = 0, in fi, 

[v £ ■ V)v E = -V/i £ , infi, 
w e • n = Vq - nln i, on dfl, 
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where n is the outward normal direction. Furthermore, as e — > 0, the corresponding vorticity 
u> £ := curlt; e satisfying 

supp (w+) C Uf =l B{z+ £ , Ce) forz+ £ G O, i = 1, • • • , m, 

supp{u-) C U] =1 B( Z r e , Ce) forzl £ G Q, 3 = 1, 



J 

v^l.e! ) Z m,Ei 1 ^n.e) ~~ ^ (^1 ! ' * * ) 4) ^1 ) * * ' ; )■ 



Proof. Since the arguments are similar to those used in section 3 and section 4 we will not 
give detail here. Let u$ be the map obtained in the reduction procedure. Define 

K(Z) = J(V+ z -V s ; z +u; s ). 

Then, as in Lemma 14.11 we can prove 

^(^C-y + o'" 1 " 11 "' 1 



\ne\P+ 2 



Similar to Theorem II A\ we can obtain a solution u £ for (15.11) . 
Define for x G f2, 

curl(w e — q £ ) 

v £ = V I H£ \ 

b | 

K u e - qs) p + L . b(-u £ - q £ y + 



h = - v . : r: + 



(p+l)e 2 (p+l)£ 2 2 

Thus, v £ is a stationary solution of ( 11. ip with 

b b 
curlv e = — (u £ - q £ ) p + zi-Ue ~ &)+• 
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Now, as in Theorem 11.11 we have 



b , _ f b 
lne\ p 



curlv e = / — [u £ -q £ f + - \ —{-u £ -q z )\ 
n Jn £ Jo, e 



[ b(v+ z -v s - z +uj s - q y 



i=l 



In 


2 


hie 




lne 





In 


2 


lne| 




lne 





j=1 * Jb (,-) v 

= y> 2?r g£ I lne | _ |l n£ | Q 

~~ 4f 6 + In 4- 6" In -4- 

-A 2^(2+) ^ 27rg(^7) 

So, the result follows. □ 

Remark 5.2. For any groen strictly local minimum points zf, ■ ■ ■ , 5~ o/ ^ on 

i/ie boundary dQ, we can also obtain the corresponding results as in Theorem ] 1. Oft 

6 Technical Estimates 

In this section we will give precise expansions of I (V$ t z), which has been used in section 4. 
Let 

^7/>*' 2 -^/^-« +1 - 

Recall that 

m 
3 = 1 

Lemma 6.1. There is a large constant L > such that 

V 5)Z (x) - q(x) < 0, x e n\ U™ l B Ls5 .( Zj ). 
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Proof. The proof is similar to Lemma A.l in [8]. For reader's convenience, we give a sketch 
here. 

If a > is small and \x — zj\ > s$ •, j = 1, ■ ■ ■ ,m, 

V s ,z - q{x) = V Sjl .j fdtg . - j^fcsfo Zj) - q{x) 

i=i V s «,i / 

, y ^ '" t g 

3=1 
m 

<5])g^t7(l + o(l))-c 

3=1 

< 0. 



If £s<5j < |x — %| < sfj, then it follows from ( 12 .11)) that 

V s ,z ~ q{x) = V 5kAiuZi - + O (jj^d^, z 



In 2 






lne 


2 



q s ,ilnL | /'In 2 | lne^ 



< 



ln^ 



/ In 2 | lng| \ 
V |lne| 2 J 



Proposition 6.2. We have 



/w ^6Sty( 1+ ^) +0 (^) 



Proof. Taking advantage of (12.61) . we find that 

\\vv 5 . z \ 2 

b 



m „ m „ 

E / / to 

3=i ^ jy< Jn 

+ IH(H) |VW+ S^( 



S,bj,q s j,Zj 



qs,j)+Vs,z,i 



□ 



24 



First, we estimate 



p+l 



bjq s , : 



S S ,j JU 

2p 2(p+l) 
2p_ / X \ — f . 2(P+1) / X \ "FT 

=^4V" — / ^ + &i4A- - 1 - . 

_ Jpl_5^ b p-i fj_\ I ( f ) P(x)g(z j + s Stj x, Zj) 

in— \s s ,jj J Bl (o) 

_ 2ix5 2 qj 3 {p+l^ qjj 27rS 2 g{z v z j ) ql j Q f ggyf^Mj 

&i 2 l ln iftl 2 *i l ln ^l 2 % V liner 1 

Next, for j ^ i, 



bj J n (V 5;ijM j>Zj -q 5j y + V s , z 
_ bjqs,i 



1,1 ^ •''<' ,,- , • 



In— ' \ss,j/ j Bl (o) 

2n5 2 qs,iqs, j G(z i ,z j ) ( \VG(z h %)l a ij 



+ 



Note that (I2.12p and || — i| = — Zj\), we can obtain that 

*7 n (H) |vw 



=0 



In d 



Similarly, 

*/ (^ WwWw = (l4) 
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By Lemma I6.1[ we have 



f b(V s , z - q y + +1 = [ b (v^Jx) - g s ,i + o ( 
Jn =1 Jb l ( Zj ) \ V 



o' p 



In 2 






hie 


2 



p+1 



In e 



Thus, we find 

m r-9 ~2 m r9 / \ "2 m r9 * 

w z) = y _ + y ™* jw* G(zj, zj) + o , 

V , ; ^-fln^-ft, ln-^ 2 6, ^^In — ln-^ ' V M 

The result follows from bj = b(zj) and the fact that 

q Sti = q{zi) ( 1 + + O ^|^| ) | , i = 1, • " • , m. 



□ 



References 

[1] G.K. Batchelor, An introduction to Fluid Dynamics, Cambridge University Press, (1967). 

[2] M.S. Berger and L.E. Fraenkel, Nonlinear desingularization in certain free-boundary 
problems, Comm. Math. Phys., 77(1980), 149-172. 

[3] R. Camassa, D.D. Holm and CD. Levermore, Long-time effects of bottom topography 
in shallow water, Phys. D, 98(1996), 258-286. 

[4] R. Camassa, D.D. Holm and CD. Levermore, Long-time shallow-water equations with a 
varying bottom, J. Fluid Mech., 349(1997), 173-189. 

[5] D. Cao, Z. Liu and J. Wei, Regularization of point vortices for the Euler equation in 
dimension two. arXiv:1208.3002, 

[6] D. Cao, Z. Liu and J. Wei, Regularization of point vortices for the Euler equation in 
dimension two, part II. arXiv:1208.5540, 

[7] D. Cao, S. Peng and S. Yan, Multiplicity of solutions for the plasma problem in two 
dimensions, Adv. Math., 225(2010), 2741-2785. 

[8] E.N. Dancer and S. Yan, The Lazer-McKenna conjecture and a free boundary problem 
in two dimensions, J. London Math. Soc., 78(2008), 639-662. 

[9] S. De Valeriola and J. Van Schaftingen, Desingularization of vortex rings and shallow 
water vortices by semilinear elliptic problem. arXiv:1209.3988, 



26 



[10] L.E. Fraenkel and M.S. Berger, A global theory of steady vortex rings in an ideal fluid, 
Acta Math., 132(1974), 13-51. 

[11] F. Flucher and J. Wei, Asymptotic shape and location of small cores in elliptic free- 
boundary problems, Math. Z., 228(1998), 638-703. 

[12] G. Li, S. Yan and J. Yang, An elliptic problem related to planar vortex pairs, SIAM J. 
Math. Anal, 36(2005), 1444-1460. 

[13] Y. Li and S. Peng, Multiple solutions for an elliptic problem related to vortex pairs, J. 
Diff. Equal, 250(2011), 3448-3472. 

[14] C.C. Lin, On the motion of vortices in two dimension - I. Existence of the Kirchhoff- 
Routh function, Proc. Natl. Acad. Sci. USA, 27(1941), 570-575. 

[15] G. Richardson, Vortex motion in shallow water with varying bottom topography and 
zero Froude number, J. Fluid Mech., 411(2000), 351-374. 

[16] D. Smets and J. Van Schaftingen, Desingulariation of vortices for the Euler equation, 
Arch. Rational Mech. Anal, 198(2010), 869-925. 

[17] J. Wei, D. Ye and F. Zhou, Bubbling solutions for an anisotropic Emden-Fowler equa- 
tion, C. R. Acad. Sci. Pairs, Sect. I, 343(2006), 253-258. 

[18] J. Wei, D. Ye and F. Zhou, Bubbling solutions for an anisotropic Emden-Fowler equa- 
tion, Calc. Var. Partial Differential Equations, 28(2007), 217-247. 

[19] J. Wei, D. Ye and F. Zhou, Analysis of boundary bubbling solutions for an anisotropic 
Emden-Fowler equation, Ann. Inst. Poincare Anal. Non Lineaire, 25(2008), 425-447. 

[20] J. Yang, Existence and asymptotic behavior in planar vortex theory, Math. Models 
Methods Appl. Sci., 1(1991), 461-475. 

[21] J. Yang, Global vortex rings and asymptotic behaviour, Nonlinear Anal, 25(1995), 
531-546. 



27 



